Non-equilibrium transport with self-consistent renormalised contacts for a 
single-molecule nanodevice with electron-vibron interaction 
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We present an application of a new formalism to treat the quantum transport properties of fully 
interacting nanoscale junctions [Phys. Rev. B 84, 235428 (2011)]. We consider a model single- 
molecule nanojunction in the presence of two kinds of electron-vibron interactions. In terms of 
electron density matrix, one interaction is diagonal in the central region and the second is off- 
diagonal in between the central region and the left electrode. We use a non-equilibrium Green's 
function technique to calculate the system's properties in a self-consistent manner. The interaction 
self-energies are calculated at the Hartree-Fock level in the central region and at the Hartree level 
for the crossing interaction. Our calculations are performed for different transport regimes ranging 
from the far off-resonance to the quasi-resonant regime, and for a wide range of parameters. They 
show that a non-equilibrium (i.e. bias dependent) static (i.e. energy independent) renormalisation is 
obtained for the nominal hopping matrix element between the left electrode and the central region. 
Such a renormalisation is highly non-linear and non-monotonic with the applied bias, however 
it always lead to a reduction of the current, and also affects the resonances in the conductance. 
Furthermore, we show that the relationship between the non-equilibrium charge susceptibility and 
dynamical conductance still holds even in the presence of crossing interaction. 
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I. INTRO 

The theory of quantum transport in nano-scale devices 
has evolved rapidly over the past decade, as advances in 
experimental techniques have made it possible to probe 
transport properties (at different temperatures) down 
to the singlc-moleculc scale. Furthermore simultaneous 
measurement of charge and heat transport through sin- 
gle molecules is now also possible^. The development of 
accurate theoretical methods for the description of quan- 
tum transport at the single-molecule level is essential for 
continued progress in a number of including molec- 
ular electronics, spintronics, and thermoelectrics. 

One of the longstanding problems of quantum charge 
transport is the establishment of a theoretical framework 
which allows for quantitatively accurate predictions of 
conductance from first principles. The need for meth- 
ods going beyond the standard approach based on den- 
sity functional theory combined with Landauer-like elas- 
tic scattering2r22 has been clear for a number of years. 
It is only recently that more advanced methods to treat 
electronic interaction have appeared, for example those 
based on the many-body GW approximatior^"— . Alter- 
native frameworks to deal with the steady-state or time- 
dependent transport are given by many-body perturba- 
tion theory based on the non-equilibrium (NE) Green's 
function (GF) formalism: in these approaches, the inter- 
actions and (initial) correlations are taken into account 
by using conserving approximations for the many-body 
self-energ}4i^"— . 

Other kinds of interactions, e.g. electron-vibron cou- 
pling, also play an important role in single-molecule 
quantum transport. Inelastic tunneling spectroscopy 



constitutes an important basis for spectroscopy of molec- 
ular junctions, yielding insight into the vibrational modes 
and ultimately the atomic structure of the junction^^. 
There have been many theoretical investigations focusing 
on the effect of electron-vibron coupling in molecular and 
atomic scale wires^i^— . In all these studies, the interac- 
tions have always been considered to be present in the 
central region (i.e. the molecule) only, and the latter is 
connected to two non-interacting terminals. Interactions 
are also assumed not to cross at the contracts between 
the central region and the leads. When electronic in- 
teractions are present throughout the system, as within 
density-functional theory calculations, they are treated 
at the mean-field level and do not allow for any inelas- 
tic scattering events. However, there are good reasons to 
believe that such approximations are only valid in a very 
limited number of practical cases. The interactions, in 
principle, exist throughout the entire system. 

In a recent paper we derived a general expression 
for the current in nano-scale junctions with interaction 
present everywhere in the systemic. With such a formal- 
ism, we can calculate the transport properties in those 
systems where the interaction is present everywhere. The 
importance of extended interaction in nano-scale devices 
has also been addressed, for electron-electron interaction, 
in recently developed approaches such as Refs. [l3ll76l |. 

In the present paper, wc also consider interactions ex- 
isting beyond the central region. We apply our recently 
developed formalism^^ for fully interacting systems to a 
specific model of a single-molecule nanojunction. Wc fo- 
cuss on a model system in the presence of electron-vibron 
interaction within the molecule and between the molecule 
and one of the leads. We show how the interaction cross- 
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ing at one interface of the molecular nanojunctions af- 
fects the transport properties by renormalising the cou- 
pling at the interface in a bias-dependent manner. We 
also study the relationship between the non-equilibrium 
charge susceptibilit}*^! and the dynamical conductance 
for the present model of interaction crossing at the con- 
tacts. 

The paper is organised as follows: In Sec. |TT1 we briefly 
recall the main result of our current expression for fully 
interacting systems. In Sec. IIII[ we present the model 
Hamiltonian for the system which include two kinds of 
electron- vibron interaction, an Holstein-like Hamiltonian 
combined with a Su-Schrieffer-Heeger-like Hamiltonian. 
In this section, we also describe how the correspond- 
ing self-energies are calculated and the implications of 
such approximations on the current expression at the 
left and right interfaces. In Sec. IIV[ we show that 
our approximations are fully consistent with the con- 
straint of current conservation. Then the effects of the 
static non-equilibrium (i.e. energy-independent but bias- 
dependent) renormalisation of the coupling at the contact 
on both the current and the dynamical conductance arc 
studied for a wide range of parameters. We also show 
that the NE charge susceptibility is still related to the 
dynamical conductance even in the presence of crossing 
interaction at the contact. We finally conclude and dis- 
cuss extensions of the present work in Sec. [V] 
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By definition Elc(w) = Vlc + (similarly for the 

CL components) where Vlc/cl S'^'c the nominal coupling 
matrix elements between the L and C regions. g£(w) 
are the GF of the region L renormalised by the interac- 
tion inside that region, where a; = r, a, < stands for the 
retarded, advanced and lesser GF components. For ex- 
ample, for the advanded and retarded components, we 

have {gt{u:))-^ = {gti^))-' " where all 

quantities are defined only in the subspace L. 

The first line in the current equation Eq. ([Ij corre- 
sponds to a generalisation of the Mcir and Wingreen2^ 
result to the cases for which the interactions are present 
in the three L,C,R regions as well as in between the 
L/C and C/R regions. The second trace in Eq. ^ 
corresponds to inelastic events induced by the inter- 
action in the L lead. When a local detailed balance 
equation helds, this terms vanishes since locally one has 

Eq. ([IJ bears some resemblance to the expression de- 
rived by Meir and Wingreen2^ when written as: 



II. GENERAL THEORY FOR QUANTUM 
TRANSPORT 



We consider a two- terminal device, made of three re- 
gions left-central-right, in the steady-state regime. In 
such a device, labelled L — C — R, the interaction — which 
we specifically leave undefined (e.g. electron-electron 
or electron-phonon) — is assumed to be well described in 
terms of the single-particle self-energy E'^^ and spreads 
over the entire system. 

We use a compact notation for the Green's function 
C and the self-energy E matrix elements M{uj). They 
are annotated Mq {Ml or Mn) for the elements in the 
central region C (left L, right R region respectively), and 
Mlc (or Mcl) and Mrc (or Mcr) for the elements 
between region C and region L or R. There are no direct 
interactions between the two electrodes, i.e. E^^^^^ = 0. 

In Refs. we showed that for a finite applied bias 

V the steady-state current Il{V) flowing through the left 
LC interface is given by: 
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where we use the definitions E^ — 
VcL 5l Vlc - i/lFl and E^'° - E^''' = Vcdgl - 
9lWlc = iFl- Hence if"^ becomes 

Hence one can see by comparing Eq. ([1]) and Eq. ([3]) 

that the quantities T lc (^lc) ^'^'^ '^lc ^-"^^ playing the 
role of the L lead self-energy E^ (S£) and E^ respec- 
tively when the interactions cross at the LC interface. In 
Meir and Wingreen model, the leads are non-interacting, 
hence the second trace Tr^^} [...] in Eq. ([T]) does not exist. 



III. MODEL FOR THE INTERACTION 

A. Hamiltonians 

We consider a single-molecule junction in the presence 
of electron-vibron interaction inside the central region 
and crossing at the contacts. Using a model system to re- 
duce these calculations to a tractable size, we concentrate 
on a single molecular level coupled to a single vibrational 
mode. A full description of our methodology, for the in- 
teraction inside the region G, is provided in Refs. [78l - l80| . 
Furthermore, we consider that the electron-vibron inter- 
action exist also at one contact (the left L electrode for 
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instance). This model typically corresponds to an exper- 
iment for a molecule chemisorbed onto a surface (the left 
electrode) with a tunneling barrier to the right R lead. 

In the following model, we consider two kinds of 
electron-vibron coupling: a local coupling in the sense 
of an Holstein-like coupling of the electron charge den- 
sity with a internal degree of freedom of vibration inside 
the central region, and an off-diagonal coupling in the 
sense of a Su-Schrieffer-Heeger-like coupling^SiSi to an- 
other vibration mode involving the hopping of an electron 
between the central C region and the L electrode. 

The Hamiltonian for the region C is 

He = £oct^d + hujQa^a + ^0(0' + a)d^d, (4) 

where (d) creates (annihilates) an electron in the 
molecular level Sq. The electron charge density in the 
molecular level is coupled to the vibration mode of energy 
ujQ via the coupling constant 70, and (a) creates (anni- 
hilates) a vibration quantum in the vibron mode luq. The 
central region C is nominally connected to two (left and 
right) one-dimensional tight-binding chains via the hop- 
ping integral tpL and ipfl- The corresponding electrode 
a ^ L,R self-energy is EJ^(a;) = tQ„//3a exp'''"''^) with 
the dispersion relation uj = Ea + 2/3aCos{ka{uj)) where 
Ea and Pa are the tight-binding on-site and off-diagonal 
elements of the electrode chains. 

To describe the electron-vibron interaction existing at 
the left contact, we consider that the hopping integral 
<0L is actually dependent on some generalised coordinate 
X. The latter represents either the displacement of the 
centre-of-mass of the molecule or of some chemical group 
at the end of the molecule link to the L electrode. At 
the lowest order, the matrix element can be linearised as 
toL{X) = toL + t'oL-^- Hence the hopping of an electron 
from the C region to the L region (and vice versa) is cou- 
pled to a vibration mode (of energy uja) via the coupling 
constant 7^ (itself related to ^ql)- The corresponding 
Hamiltonian is given by 

Hlc - lA{b^ + b){cld + Scl) + uJAb^b, (5) 

where (&) creates (annihilates) a vibration quantum in 
the vibron mode u:a, the generalised coordinate is X = 
^/K/J2rrlA0OA){b'' + b), and (cl) creates (annihilates) 
an electron in the level El of the L electrode. 

The Hamiltonians Eq. (U) and Eq. (O are used to cal- 
culate the corresponding electron self-energies at differ- 
ent orders of the interaction 70 and "/a using conventional 
non-equilibrium diagrammatics technique o^^'^^ . 

Furthermore, at equilibrium, the whole system has a 
single and well-defined Fermi level 11°^. A finite bias 
V, applied across the junction, lifts the Fermi levels as 
= + ''lL,ReV. The fraction of potential drop^ 
at the left contact is 77^ and ?7i? = 77^ — 1 at the right 
contact, with /i^ — fi^ = eV and rji G [0, 1]. 



B. Self-energies for the interactions 

The electron-vibron self-energies in the central region 
C are calculated within the Born approximation. The de- 
tails of the calculations are reported elsewhere^i^ so we 
briefly recall the different expressions for the self-energies 

U = + Zjp With 

E^- = S^- = 2^ / ^^G<iu;') = -2:^ {no) , (6) 
ujf) J Zn Wo 

with (nc) = ~i J diu/2Tr Gq{uj) and 

I]^'^(w) = 77o' I ^ Dfiu) G%{u u) , (7) 

and 

E^''-(w) = J ^Dl{u - u) (G<(u) + G^u)) 
+D<iuj^u)G},iu) , 
with the usual definitions for the bare vibron GF Dq-. 

Df{u) = -2m [{No)S{uj T ^^o) + (No) + l)5(c^ ± c^o)] 



(9) 

where (Nq) is the averaged number of excitations in 
the vibration mode of frequency ujq given by the Bose- 
Einstein distribution at temperature T^ih- In the follow- 
ing, we work in the limit of low temperature for which 
{No) = 0. 

As a first application of our transport formalism for 
crossing interactions, we consider a mean-field approxi- 
mation for the electron-vibron coupling at the LC inter- 
face. This leads to the Hartree-like expressions for the 
many-body self-energies at the LC interface: 

S^c'''/" = -2^(n,.c), (10) 

UJA 

where 

Kc) = ^ Gtc(^)- (11) 

Similarly the self-energy E^]^''^^" is obtained from 
("cl) = / da;/27r G^^(w). 

One can see that the interaction crossing at the LC in- 
terface induces a static (however bias-dependent) renor- 
malisation of the nominal coupling Vol = Vlc = toL 
between the L and C regions. This non-equilibrium 
renormalisation will induce, amongst other effects, a bias- 
dependent modification of the broadening of the spectral 
features of the C region. Since the renormalisation is 
static at the Hartree level, we use a "potential" notation 
to represent the normalised coupling: E^^^ = Vlc with 

Vlc = VLc-2^{nLc), (12) 

UJA 
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and similarly for Vcl- ^ 0.1 to ~ 5.0. One should note, however, that in all 

The static renormalisation of the nominal coupling tpL the calculations we have performed, the density matrix 

is driven by the ratio Jj^/lua- In the following numerical element {ulc) is always of the order of 10"^. Therefore 

applications, we consider small to larger renormalisation the renormalisation effects are always smaller than the 

effects, for which the ratio 7^/(w^ioL) is ranging from nominal coupling ioL itself. 



we need the off-diagonal elements and G^^. 



expression for the G^^ GF matrix element is obtained from the corresponding Dyson equation G^^ = [gSGJ^^. 



In order to get the renormalised couplings Vlc/cLj 

nent if 

After formal manipulation, we find that 

G<c{u) = g<VLcGi + Glc {VcLgtVLC + + ^^ci?5^^i?c) G^ 
with Gp is the renormalised advanced GF of the central region 

Gl[u) = [[g'^cV -^i- VcRgRVRc - VcLglVLc 
and the off-diagonal element G'^^(j is given by 

Glcioj) = GIVlc [WcV - VcBg^RVncY' , 



The closed 



(13) 



(14) 



(15) 



where g^ is the renormalised GF of the disconnected C region [p^] ^ 
we get the following expression for G^^ by using the notation Y^'^ - 
self- energy: 



: [g^] ^ — Sp. After further manipulation, 
ycRgli^RC for the non-interacting R lead 



G 



LC 



[9l\ 



Vlc 



Y, 



R,r 
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^VcL 
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¥c]-' 



Y, 



(16) 



Similarily we can calculate the off-diagonal element G^^ from the corresponding Dyson equations G^^ — [gEGJ^^. 
We find 



G<^iuj) = G^cVcLgt + G^c [VcLglVLc + S< VcRg^Vnc Gl 



where Gp(w) is the retarded version of Eq. and 

G"cl(^) = \Wc 



Y, 



R,a 



'^VclGI. 



^CLi 



(17) 



(18) 



As expected from the definition of the different GFs, we can see that indeed (Gp^(a;))* = — G^p(aj) and 
(G^Jc^))* = Glciuj). ^ 



C. Calculations 

Calculations are performed in a self-consistent manner. 
There are different ways to calculate the GF and the self- 
energies in a self-consistent way: in the present work, we 
first perform self-consistent calculations for the central 
region G. The new hopping matrix elements between the 
G and L regions given by Eq. are then updated via 
the use of Eqs. ([T7)) and p^ . The calculations are then 
re-iterated until full self-consistency for all self-energies 
and GF in the region G and at the left contact E^^J^^^^ 
is obtained. 

For the model of crossing interaction we are consider- 
ing here, further simplications can be introduced in the 
calculation of the current. Since there is no other inter- 



action inside the L and R regions, the current expression 
is given only by the first line of Eq. ((!]). Furthermore, 
there are no lesser and greater components for the self- 
energy at the mean-field level. The current expres- 



'LC 

sion Eq. jH) 



reduces to: 



e f du 



Tr 



{C} 



g: 



c'^c 



(t^)t) 



(19) 



with the following simplified expressions for the Tc quan- 
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titles 



^ c 
(t^)t 



C 



VcL gl Vlc, 
VcL gl Vlc, 

VcL gt ^LC, 

-VcL gt Vlc, 



(20) 



since E^'^ = Vlc and = 0. 

Furthermore, the current at the CR interface has a 
conventional Meir and Wingreen expression, since there 
are no crossing interactions: 



rMW 



2^ 



Tr 



{C} 



R.< 



^c 



(21) 



where we use the conventional definition for the R lead 
self-energies: Ti^''^ {lo) = Y^'''{uj) = Vc Rg%{i^)Viic with 
VcR = Vrc = toR. 

The main differences between the expression for the 
left and right currents and arise from the non- 

equilibrium static renormalisation of the coupling be- 
tween the central region C and the L lead given in 

Eq. dni). 



IV. RESULTS 

Wc have perfomed calculations for many different val- 
ues of the parameters in the Hamiltonians. We present 
below the most characteristic results for different trans- 
port regimes and for different coupling strengths 7a, 
while the interaction in the region C is taken to be in the 
intermediate coupling regime 70 /wq = 0.7. The nominal 
couplings between the central region and the electrodes 
toL,R, before NE renormalisation, are not too large, so 
that we can discrimate clearly between the different vi- 
bron side-band peaks in the spectral functions. The val- 
ues chosen for the parameters are typical values for real- 
istic molecular junction o^°i^^ . 

The different transport regimes considered bellow are 
called the far off- resonant, the off- resonant, the interme- 
diate regime and the quasi-resonant regime. They corre- 
spond to different position of the molecular level Eq with 
respect to the Fermi level at equilibrium. With the 
static renormaliscd molecular level Eq — £q — Jq/ujo, the 
far off-resonant regime corresponding to Sq — fj,'^'^ ^ 
and ^ several wq, the off- resonant regime corresponds 
to £0 ~ /^'^'^ S> and ^ one or two wo, the intermedi- 
ate regime corresponds to eo — ^ ^0 and the quasi- 
resonant regime corresponds to Eq — /i""^ < ljq. 

In the following the current is given in units of charge 
per time, the conductance in unit of quantum of conduc- 
tance Go = 2e'^/H and the bias V and the normalised 
coupling VcL in natural units of energy where e = 1 and 
h=l. 




FIG. 1: (Color online) Current at the left and right interfaces 
Il{V) (black line with dots) and Ir{V) (red line) given re- 
spectively by Eq. H19p and Eq. (|21|) . The current is fully con- 
served as can be seen by the green dashed line which shows 
II + Ir magnified by a factor of 2000. The parameters are 
£0 = 0.3, ujo = 0.30, 70 = 0.21, toR = toL = 0.15, 7a = 0.11, 
LJA = 0.10, l3a = 2.0,60 = 0.0. 



A. Current conservation 

One of the most important physical conditions that our 
formalism needs to fulfil is the constraint of current con- 
servation. Wc use conserving approximations to calculate 
the interaction self-energies in the central region G and 
for the crossing interaction at the left interface. However, 
since there is no interaction crossing at the CR inter- 
face while the interaction is crossing at the LC interface, 
we have to check that the current given by Eq. ([T^ for 
/l(I^) is equal to the Meir and Wingreen current given by 
Eq. dill) for Ir{V), i.e. Il + Ir = 0. Figurcfflshows that 
the condition of current conservation is indeed fulfilled, as 
expected. We have carefully checked that the current is 
conserved for all the calculations presented in the present 
paper, i.e. that \Il + Ir\/\Il\ - |/l + < 10"^ 



B. Static non-equilibrium renormalisation 

In Figures [H to [51 we show the dependence of the cur- 
rent I{V), of the dynamical conductance G{V) = dl/dV 
and of the renormaliscd coupling Vcl{V) on the applied 
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FIG. 2: (Color online) Far off-resonant transport regime with 
£o ~ 0.9. Current I{V) (left panel), dynamical conduc- 
tance G{V) = dl/dV in logarithmic scale (top-right panel) 
and the real part of non-equilibrium renormalized static cou- 
pling Vcl{V) (bottom-right panel) for different values of 
the coupling 7a strength at the left interface between 0.0 
and 0.14, corresponding to the follow ratio: 'y'^/(ujAtoL) ~ 
0.0,0.167,0.427,0.807,0.96,1.127,4.667. The NE renormal- 
ization of Vcl{V) corresponds to a non-monotonic reduction 
of the nominal hopping integral toL with applied bias, which 
induces a diminution of the current, and modifications of the 
widths of the conductance peaks accompanied with a slight 
shift of the conductance peaks at high bias. See text for a 
complete discussion. The other parameters are ujq = 0.30, 
70 = 0.21, toR = toL = 0.15, 7A = 0.10, /3„ = 2.0, = 0.0. 



bias V, for different values of the interaction strength 
7A at the left contact. Wc consider different transport 
regimes ranging from the far off- resonant regime (Fig. [5]) , 
the off- resonant regime (Fig.|31), the intermediate regime 
(Fig. m and finally the quasi-resonant regime (Fig. [5]). 

The NE renormalisation of the coupling at the left con- 
tact given by Eq. (fT2l) corresponds to an effective de- 
crease of the hopping integral leading to a decrease of 
the current for increasing values of ja, as can be seen in 
the left panels of figs. [2]to[S] The real part of Vcl{V) 
shows a non-monotonic behaviour with the applied bias 
and presents features (local dips) at applied biases corre- 
sponding to peaks in the dynamical conductance (bottom 
right panels of figs. [5] to [5]). Therefore, the interaction 
crossing at the left contact not only decreases the value 
of the current but also affects the width of the peaks in 



FIG. 3: (Color online) Intermediate off-resonant transport 
regime with eo = 0.5. The other parameters are the same as 
in figure [21 See text for a complete discussion. 



the conductance, as can be seen in the top right panels 
of figs. H to El 

In all our calculations, it appears that each conduc- 
tance peak now has an asymmetric shape, i.e. a different 
broadening on each side of the peak, which is due to the 
non-monotonic and asymmetric behaviour of Vc l versus 
applied bias. The detailed understanding of such a be- 
haviour is rather complex. However, one can obtain a 
qualitative understanding of our results by considering 
the following analysis. 

Wc can use the current expressions for 1^ and Ir given 
by Eq. (|T^ and Eq. respectively to obtain a sym- 
metrized current / = [II — Ir)/'^ expression within a 
good approximation^^ : 



where 7rAc(w) = — 3mG'^(w) and 



(22) 



(23) 



with the self-energy in the central region 1^^^^'^ (lu) given 
by the sum of Eq. ^ and Eq. ^ ; and we recall that 



Sp''(a;) = VcL gl{i^) Vlc 
VcR 9r{^) Vrc ■ 



R.r 



zr«/2. 



(24) 
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FIG. 4: (Color online) Intermediate off-resonant /resonant 
regime with eo = 0.3. The other parameters are the same 
as in figure O See text for a complete discussion. 



FIG. 5: (Color online) Quasi-resonant transport regime with 
£o = 0.2. Same parameters as in figure [J] See text for a 
complete discussion. 



In the off-resonant regime and for small bias V < uio, 
one can consider that 3mS^^^''^ ^ (see Ref. [73l|). 
hence 



did 



2tt {u 



rg)V4 



(25) 



We can now use Eq. (^5]) to understand the behaviour 
of the current. In Figure|6]we show the parametric curves 
I{V) — Vc'l{V) obtained for the far off-resonant transport 
regime (shown in Fig. [5]) . Figure IHl shows a complex de- 
pendence of the current versus the NE renormaliscd cou- 
pling VcL- However at low bias, one can consider that 
Fp and F^ are independent of oj and hence take such 
quantities out of the integral in Eq. (j25[) . Therefore we 



have / cx F^ and the current is quadratic in Vcl since 
Fp cx [Vchf' ■ Such a dependence can be clearly seen in 
Figure ini for the low bias regime where II < 0.01. 

At larger bias, one can no longer neglect the uj depen- 
dence of F^ (and F^) and more importantly the effects of 
the interaction in the central region, i.e. 
Hence the quadratic dependence of / on Vcl is lost. 

We now turn onto the effect of the strength of the 
nominal hopping integral toL on the renormalised cou- 
pling Vcl{V) at fixed values of the crossing interaction 
strength 7^. Figure [7] shows the relative dependence of 
Vcl{V) on the nominal coupling ioL at the left inter- 
face versus applied bias. The figure shows that there is 



no simple relationship between VchiV) and Iql for the 
whole range of parameters explored. There is a small lin- 
ear regime at low applied bias, otherwise the dependence 
of Vcl on both toa and V is highly non-linear. There is 
a progressive washing-out of the features in Vcl{V) for 
increasing values of toaj since a general increase of the 
coupling to the leads generates to a global broadening of 
the features in the spectral functions. 



C. Non-equilibrium charge susceptibility 

In a recent paper—, we have developed the concept 
of the generalized susceptibilities for nonlinear systems24 
and applied it to the charge transport properties in two- 
terminal nano-devices. We have introduced the non- 
equilibrium charge susceptibility x?^(^) concept, de- 
fined by x^^{V) = d{nc) / dV . We have shown that 
X^^(l^) is related to the dynamical conductance G{V). 
The relationship is formally different than the one ob- 
tained at equilibrium. In spectroscopic terms, both 
^^{V) and G{V) contain features versus applied bias 
when charge fiuctuation occurs in the corresponding elec- 
tronic resonance. This relationship has been demon- 
strated for model calculations of interacting nanoscale 
devices but only when the interaction is present in the 
central regional. We now check the validity of such a re- 
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FIG. 6: (Color online) Parametric curve I{V) — Vcl{V) for the 
far off-resonant transport regime shown in Fig. [2] (using the 
same legend for the different 7a values). The starting points 
of the parametric curves, i.e. V — 0, axe at the bottom of the 
graph (zero current values), and increasing V corresponds to 
increasing values of the current. At low bias {II < 0.01), II 
is quadratic in Vcl- 




V Bias 

FIG. 7: (Color online) Renormalised coupling Vcl{V) for the 
off-resonant transport regime (eo = 0.50) and 7.4 = 0.14 ver- 
sus applied bias and for different values of the nominal hop- 
ping integrals toa (a = L,R). The relative dependence of 
Vcl, shown here for (KeVci(T^) — toa)/toa, is linear only for 
small bias. Otherwise the dependence of Vcl on both toa and 
V is highly non-linear. The other parameters are the same as 
in Fig. H 



lationship between x?^(^) ^-^d G{V) when interaction 
crossing at the left contact is also taken into account. 

Figure [5] shows the results we obtain for X^^i^) ^^'^ 
G{V) for different sets of parameters corresponding to 
the ofF-resonant and quasi-resonant transport regimes. 
Once more, we find that x^^(^) ^^'^ C!{V) present 
features at the same applied bias. Such a behaviour 
validates the relationship between the non-equilibrium 
charge susceptibility and the dynamical conductance 
even when interactions cross at the contacts (at least at 
the static mean-field level). 

Finally, one should note the negative contribution to 
X^^iV) in panel (c) of Figure El Such a behaviour orig- 
inates from the properties of ^^nd not from the ap- 
proximation used to calculate the crossing interaction!^^. 
In fact, there are always two contributions to x^^, one is 
positive and corresponds to electron fluctuation and the 
other is negative and corresponds to hole fluctuation. By 
electron (hole) fluctuations, we mean the variation of the 
occupancy (versus applied bias) of electronic resonances 
located nominally above (below) the Fermi level at equi- 
librium (or at small bias). Therefore, for the off-resonant 
transport regime with Eq ^ ii'^'^, the features in x^^(^) 



corresponds to positive peaks (as shown in Panels (a ) an d 
(b) in Figure [5] and in the different figures of Ref. [Tg ) ■ 
For the off-resonant regime with Eq <^ /x®'^, one would 
get negative peaks in Xc^i^)- For the intermediate and 
quasi- resonant transport regime, one obtains both posi- 
tive and negative contributions in X?^(^) shown in 
panel (c) of Figure [51 The most peculiar case corre- 
sponds to a fully electron-hole symmetric system, i.e. 
the spectral function Ac{uj) is fully symmetric around 
the equilibrium Fermi level and for any applied bias V 
with idL.B^ = M°'' i eV/2. In that case, one can easily 
show that (tic) = —i J diu Gq{u!)/2tt is actually given 
by {nc) ^ I-x '^'^ ^c(w) for symmetry reasons. Hence 
(tt-c) is independent of the applied bias for conserving 
approximations for the self-energies, and consequently 
X^^(F) = 0. Such a behaviour can also be interpreted 
with the previous picture: for a fully electron-hole sym- 
metric system, any contributions from electron fluctua- 
tion is exactly cancelled out by the opposite contribution 
from hole fluctuation, and X?^(^) is flat and equal to 
zero for each bias^. 
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FIG. 8; (Color online) Non-equilibrium charge susceptibility 
X^^{V) and dynamical conductance G{V) versus applied bias 
for different sets of parameters. Xc'^(^) is rescaled by a fac- 
tor of 7. (a) Off-resonant regime eo ~ 0.5 and -ja ~ 0.08. 
(b) Off-resonant regime eo ~ 0.5 and ja = 0.14, (c) Quasi- 
resonant regime eo = 0.2 and ja = 0.05. x^^i"^) and G{V) 
present features at the same applied bias. Note the negative 
contribution to around V ~ 0.56. See main text for a 
complete discussion. The other parameters are the same as 
in Fig. H 



V. CONCLUSION 

We have studied the transport properties through a 
two-terminal nanoscale device with interactions present 
not only in the central region but also with interaction 
crossing at the interface between the left lead and the 
central region. To calculate the current for such an fully 
interacting system, we have used our recently developed 
quantum transport formula— based on the NEGF for- 
malism. As a first practical application, we have con- 
sidered a prototypical single-molecule nanojunction with 
electron-vibron interaction. In terms of the electron den- 
sity matrix, the interaction is diagonal in the central re- 
gion for the first vibron mode and off-diagonal between 
the central region and the left electrode for the second vi- 
bron mode. The interaction self-energies are calculated in 
a self-consistent manner using the lowest order Hartree- 
Fock-like diagram in the central region and the Hartree- 
like diagram for the crossing interaction. Our calcu- 
lations were performed for different transport regimes 
ranging from the far off-resonance to the quasi-resonant 



regime, and for a wide range of parameter values. They 
shown that, for this model, we obtain a non-equilibrium 
(i.e. bias dependent) static (i.e. energy independent) 
renormalisation of the nominal hopping matrix element 
between the left electrode and the central region. The 
renormalisation is such that the amplitude of the cur- 
rent is reduced in comparison with the current values ob- 
tained when the interaction is only present in the central 
region. Such a result could provide an partial explana- 
tion for the fact in conventional density-functional based 
calculations, the values of the current are always much 
larger than in the corresponding experiments, since no 
non-equilibrium renormalisation of the contacts is taken 
into account in those calculations. However, even though 
it provides the right trends, the decrease in the current 
obtained by NE renormalisation of the coupling to the 
leads is not as important as the effects obtained from a 
proper calculation of the band-gap and band-alignment 
in realistic molecular syste m^'^i^^" — . 

The NE static renormalisation of the contact is highly 
non-linear and non-monotonic in function of the applied 
bias, and the larger effects occur at applied bias corre- 
sponding to resonance peaks in the dynamical conduc- 
tance. The conductance is also affected by the NE renor- 
malisation of the contact, showing asymmetric broaden- 
ing around the resonance peaks and some slight displace- 
ment of the peaks at large bias in function of the coupling 
strengh ja- 

Furthermore, we have also shown that, even in 
the presence of crossing interactions, the relationship 
between the NE charge susceptibility and dynamical 
conductancoSi still holds for the different transport 
regimes considered here. 

Extensions of the present study are now considered. 
One route is to develop more accurate NE renormali- 
sation by considering, for example, a quasi-particle ap- 
proach within a dynamical mean- field- like treatement of 
the crossing interaction self-energji^. Another route is to 
go beyond the Hartree approximation for the crossing in- 
teraction self-energy using other many-body diagrams^. 
This would lead to dynamical NE renormalisation of the 
contact involving inelastic scattering processes. 

Finally, it should be noted that for the model system 
we considered here, one could also solve the problem by 
using more standard approaches (typically the original 
Meir and Wingreen approach) by extending the size of 
the central region to include all the interaction. The re- 
sults then obtained should be strictly equivalent to our 
calculations. We have already commented in detail on 
this point (in a formal and theoretical point of view) in 
Refs. [73H73 |. However our approach offers a more intu- 
itive physical interpretation of the results, i.e. renormal- 
isation of the contact (in the present case in a static NE 
mean- field scheme). Using an extended central region 
will not provide an easy physical interpretation of the re- 
sults; and as a matter of principle, it will not always be 
possible to increase at will the size of the central region, 
more especially when one considers future application of 
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the method to much large molecular system with much more complex coupling to the leads. 
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